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Abstract 

The nuclear breathing-mode giant monopole resonance is studied within an improved relativistic 
Boltzmann-Uehhng-Uhlenbeck (BUU) transport approach. As a new feature, the numerical treat- 
ment of ground state nuclei and their phase-space evolution is realized with the same semiclassical 
energy density functional. With this new method a very good stability of ground state nuclei in 
BUU simulations is achieved. This is important in extracting clear breathing-mode signals for the 
excitation energy and, in particular, for the lifetime from transport theoretical studies including 
mean- field and collisional effects. 

PACS numbers: 21.65.-f, 21.60.-n, 24.30.Cz, 21.60.Ev 
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I. INTRODUCTION 



A giant monopole resonance (GMR), i.e. a collective isoscalar O"*" excitation of a nucleus, 
has been extensively investigated in the past few decades, both theoretically and experi- 
mentally (see Ref. for the most recent review). It is nowadays well established, that 
the GMR is a nuclear compression mode governed mostly by the incompressibility modulus 
of nuclear matter ^ , which is the key quantity for the description of nuclei, supernovae 
explosions, neutron stars, and heavy-ion collisions. 

Many microscopic models have been developed for the theoretical description of giant 
resonances. They can be divided into two groups: purely quantum mechanical approaches 
and semi-classical dynamical models. 

The former group includes the models based on the nonrelativistic Isl kl or relativistic 

n n 

[5|, |6[ static Hartree-Fock plus random phase approximation (RPA) method. The RPA 
technique can be derived from a more general time-dependent Hartree-Fock (TDHF) theory 
(c.f. [71, l8|) in the case of small-amplitude excitations. The damping of a collective mode in 
a pure mean field RPA picture originates from the coupling to the particle-hole excitations 
(Landau damping or fragmentation width) and from coupling to the continuum states, which 
is equivalent to the particle loss in TDHF calculations [7|. Also constrained relativistic 
mean-field approaches have been developed and applied in the case of a GMR IJIlol. The 
collective nature of the GMR in these quantum mechanical prescriptions manifests itself 
from a coherent superposition of many single-particle transitions from one major shell to 



another 



Generally, within uncertainties of an underlying energy-density functional. 



the fully quantum approaches are able to describe the GMR energies for various sets of 
nuclei. However, the conclusions of different authors on the GMR total width are clearly 
different. Sagawa et al. [4] claim that the Landau damping and coupling to the continuum 
states explain the major part of the total GMR width for the ^osp^ nucleus and Sn isotopes. 
However, Piekarewicz and Centelles state that the RPA calculation fails to account for the 
spreading component of the full escape-plus-spreading width because of the lack of coupling 



to the more complex than partic 
The second group of models 



e-hole configurations. 



12| - 21| solves the BUU(-like) or Vlasov(-like) equations. 



which are the semiclassical limits of a quantum kinetic equation [22]. An advantage of the 
kinetic transport equation with respect to the TDHF theory is that binary collisions are 
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naturally included, apart the interaction between the particles due to the classical nuclear 
mean-field. In the classical picture, a GMR can be qualitatively understood in terms of a 
radial collective vibration of a nucleus: protons and neutrons oscillate in phase with a certain 
amplitude, which is damped due to dissipation. For this reason, the GMR is often referred 
to as a "breathing mode". The frequency of the breathing mode characterizes the excitation 
energy and the temporal damping of the amplitude — the life time of the resonance. 

Semi-classical treatments of a GMR in finite nuclei have been so far restricted to pure 



Vlasov dynamics 



18 



20 



2l|, where any collisional effects are completely neglected. Thus, 



these models predicted the excitation energy of the breathing mode appropriate well, but not 

assical description of the GMR 
17|, which, however, takes into 



the width. Alternatively, the collisional effects in the semi-c 
were modeled within a linearized Landau- Vlasov equation 
account the finite size effects only very approximately. There are several works where the 
Boltzmann-Nordheim- Vlasov approach - which takes into account the selfconsistent nucleon 
mean field and nucleon- nucleon collisions - has been applied to the collective dipole 23|, l24 1 
and quadrupole 25[ motions excited in heavy-ion collisions at low energies. 

The present work is an attempt to describe simultaneously the centroid energy and the 
width of a GMR in finite nuclei. To this aim, we perform the full BUU calculations taking 
into account both the mean field and collision term. The results of the full BUU calcula- 
tions are then compared with the results of solution of the Vlasov equation to make the 
quantitative conclusions on the contribution of two-body collisions to the total GMR width. 

The numerical solution of a BUU equation for the case of a small- amplitude collective 
vibration excited in a finite nuclear system is extremely difficult. It requires a very good 
stability of the ground state configurations, which is difficult to reach in a test-particle 
technique underlying any numerical method to solve the BUU equation. So far empirical 
density distributions have been used to initialize ground state nuclei in transport theoretical 
simulations, which might be not always consistent with the energy density functional used for 
;he propagation of the system. Another well known problem (see e.g. discussions in Refs. 
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14j ) is related to the calculation of Pauli blocking factors in the Uehling-Uhlenbeck 



collision integral for the small-amplitude Fermi surface distortions 

We thus have improved the relativistic transport approac 
Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model 



26|-|28| based on the Giessen 
29| by performing relativistic 



Thomas-Fermi (RTF) calculations with the same energy-density functional as that used 
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in the dynamical evolution. In this context, the isovector- vector p meson field together with 
its gradient terms have been included in the calculations of the present work. The initializa- 
tion of neutron and proton densities according to the RTF calculation largely improves the 
ground-state stability in numerical simulations of the Vlasov and BUU dynamics. We have 
also worked out to improve the numerical treatment of the Pauli blocking in test-particle 
simulations of the small-amplitude nuclear motions, in-particular, in the nuclear surface 
regions. 

The structure of the work is as follows: the standard theoretical background is presented 
in section [TTl The modified initialization procedure and numerical treatment of the Pauli 
blocking are then presented in section IIIII In section \YV\ results of the ground state simula- 
tions are discussed, before the calculations of the GMR excitations are presented. Finally, 
conclusions and outlook are presented in the last section |Vl 



II. THE RELATIVISTIC TRANSPORT EQUATION 



The baryonic mean-field is modeled within the non-linear Walecka model (mean-field 



approximation of the QHD) 30|,l311- The non-linear Walecka model Lagrangian includes the 



nucleon field ip, the isoscalar-scalar a meson field, isoscalar-vector uj meson field, isovector- 
vector p meson field and the electromagnetic field A and reads {h = c= 1) as: 
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A^) - (m + g,a)]i; 



(1) 



where Vt^^ = d^ujy - dyU^, R^y = d^p^ - d^p^ and F^^ = d^Au - dyA^ are the field 
tensors. In Eq. ([1]), the arrow above a symbol indicates the isovector character of the 
corresponding field. The term U (a) = ^m'^a^ + \g2C^ + jdsO'^ contains the selfinteractions 
of the a field, added according to 



31| . The bare hadron masses m, m^, and rrip, 



coupling constants g^, guj, dp and the non-linear parameters g2 and g^ have been adopted 
from the NL3* parametrization of the non-linear Walecka model, which gives reasonable 
values for the incompressibility modulus, Koo = 258 MeV, and the nucleon Dirac effective 
mass m* = 0.594m at the saturation density po = 0.150 fm~^ of nuclear matter 32 |. 



NL3* parametrization is the modification of the well known NL3 set of parameters 



he 



33| 



adjusted to describe the ground state properties of both spherical and deformed nuclei as 



well as the GMR energies of heavy nuclei, 
optical potential, studied recently in Ref. 



he momentum dependence of the proton-nucleus 



is of minor importance here. 



The theoretical description of heavy-ion collisions is realized within the GiBUU transport 



approach 29j, which is based on a relativistic kinetic equation. Thorough derivations of the 



Tansport equation from an effective hadron- meson field theory SOj, can be found elsewhere 



35| . The relativistic kinetic equation reads 

j=n,p 

with yU = 0, 1, 2, 3 and k = 1,2, 3. The l.h.s. of Eq. ([2]) describes the classical Vlasov prop- 
agation of a one-body phase space distribution function fi{x,p*) for protons and neutrons 
{i = p, n) in the mean meson fields. This is expressed in terms of a kinetic four-momentum 
p* = Pi—Vi, where pi is the canonical four-momentum, of the field tensor F/"^ = d^V^'^—d'^Vl^, 
and of the Dirac effective mass m* = m + S. Here Vl" and S are the vector^ and scalar field, 
respectively. 

Vt = g^uj^ + gpT^p'^ + I (1 + ri^)^^ (3) 
S = g.a , (4) 

where rf = +(— )1 for i = p{n). 

The particles are assumed to be on the Dirac effective mass shell, i.e. 



p " = ^/rn*^~+~p*^ . (5) 

In fact, the GiBUU model propagates not only nucleons, but also all resonances up to 
the mass of 2 GeV, as well as mesons, e.g., pions, kaons. However, in the present study 
we will concentrate only on the nucleonic degrees of freedom. In this case, the collision 
terms in the r.h.s. of Eq. (|2]) are the Uehling-Uhlenbeck collision integrals describing elastic 
nucleon-nucleon scattering: 



^ Assuming that there is no mixing between proton and neutron states, the 1-st and 2-nd isospin components 
of the p-meson fields vanish. 



5 



where the time-space argument x of the distribution functions is dropped for brevity and 
= p*. The hole distribution functions are denoted as /j(p*) = (1 — /i(p*)). Thus the 
final state Pauli blocking is explicitly included in (jS]). The collision integral (jS]) depends 
on the differential elastic scattering cross section rfcTy (p*, pg; Pi/, P2/) with initial momenta 
pj, P2 and final momenta p]^,, P2/, and on the relative velocity f 12 of colliding nucleons. We 
use in this work the energy- and angular-dependent vacuum nucleon-nucleon cross sections 



(see |26| for details). 

An exact solution of the set of the coupled transport equations for the different hadrons 
is not possible. Therefore, the commonly used test-particle method for the Vlasov part 
is applied, whereas the collision integral is modeled in a parallel-ensemble Monte-Carlo 
algorithm. The GiBUU transport model has not been applied to low energy reactions so far. 
An important requirement here is a very good description of the ground state of a nucleus 
in test particle method, which is the topic of the next chapter. 

III. IMPROVED METHOD FOR NUCLEAR GROUND STATES IN GIBUU 



Practical 
nique (c.f. 



ly all existing nuclear kinetic transport models based on the test particle tech- 



20 



21 



23 



29| and refs. therein), including GiBUU, use the same method of the 
nuclear ground state preparation. Typically, the coordinates of test particles are sampled 
according to empirical Woods-Saxon or even uniform density profiles, while the test particle 
momenta are distributed with the help of a local density approximation. The standard nu- 
merical treatment of the transport equation works well for high energy reactions, in which 



several collective flow observables are described quantitatively well |36|. However, in low 
energy reactions the memory of the exit channel to the initial configuration is important. A 
disadvantage of standard numerical treatments is that the initial distribution functions of 
protons and neutrons deviate from the corresponding static solutions of the Vlasov equation. 
Therefore, nuclei are not initialized in their proper ground states. This makes their dynam- 
ical propagation unstable. Another source of spurious instability is the numerical treatment 
of the Pauli blocking, which mainly infiuences the momentum space of the test particles in 
full BUU simulations. 
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A. Phase space initialization 



We improve the phase space initiahzation of ground state nuclei in the transport model 
in the following way: the nuclear ground state is described in a semi-classical treatment, in 
which the density distribution of a spherical ground state nuclei is obtained by minimizing 
the energy functional. In RMF, the energy functional corresponds to the relativistic Hamil- 
tonian density, which is obtained from the T°'^-component of the energy-momentum tensor. 
The same functional is then used for the propagation of the system. In calculations of the 
time evolution, we neglect the time derivatives of the mean mesonic fields and of the electro- 
magnetic field in the Lagrangian ([T]). For the mean mesonic fields, this should be a reasonable 
approximation, since the time scale of the GMR motion ~ 40 — 80 fm/c (c.f. Fig. [7] below) 
is much larger than the time scale of the free field oscillations ~ l/mmes = 0.2 — 0.4 fm/c, 
where mmes ~ 0.5 — 0.8 GeV is a meson mass. For the electromagnetic field, neglecting time 
derivatives corresponds to disregarding a radiation. The last could be, in-principle, treated 
in terms of a quantum transition probability. However, in the specific case of 0"^ transitions, 
the electromagnetic radiation is suppressed js^]- The space components of electromagnetic 
field are also neglected, i.e. A = (^°, 0,0,0). In the most calculations, in order to save 
CPU time, we also drop the space components of the vector meson fields. However, for 
completeness, the formalism below takes into account these components. 

The RMF Hamiltonian density of the non-linear Walecka model reads: 



e = r 



•00 



\3 



J2 J rfyrf(x,p*)/.(x,p*) + ^(Va)2 + f/(a) 



(27r 

' i=p,n 

- i [Vu;^ Vu;^ + m>2 + V p^^^V p\ + (p3)2 + ( v^O)^] . (7) 

The Hamiltonian density (171) takes into account the gradients of the meson fields in 
coordinate space according to |27|. This is important for description of surface effects, 
which is impossible in the usual local density approximation neglecting the gradient terms 



fc.f. 



21 



26|). In distinction to the previous GiBUU calculations in the RMF mode, the 



isovector p-meson field and the Coulomb field ^0 are explicitly included now. 
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The meson and electromagnetic field equations have the following form 

(-A + ^).^-..P., (8) 

{-/\ + ml)uj^ = g^f^ , (9) 

{-A + ml)p'^ = g,jf , (10) 

-A^'^ = ej; . (11) 

The source densities and currents in the r.h.s. of Eqs. fl8l)- f|TTl) are expressed in terms of the 
distribution functions as follows 

Psi^) = T^E /^V^/.(x,P*), (12) 

i=p,n 

2 f p*^ 

ji(^) = J dY-^fi{x,P*) , ii=P,n) (13) 

fsi^) = J^ix)+ji:ix) , (14) 

jT(^) = • (15) 

In the static case, the distribution functions are Fermi distributions: 

/r^^(r,p*) = 0(p^^(r)-|p*|). (16) 

The space components of the source currents ( fT3|) -(fT5l) as well as those of the mesonic and 
electromagnetic fields disappear in a static system. Therefore, the Hamiltonian density 
becomes a functional of proton and neutron densities pp,„ = „ only. 

The RTF equations for a static nucleus with Z protons, neutrons and A = N + Z nu- 
cleons are obtained by applying a variational principle to the total energy E = f e[pp, Pnjfi^r 
under the constraint of particle number conservation: 

^ j (e[Pp, Pn] - PpPp(r) - p„p„(r)) £r = . (17) 

The chemical potentials for protons and neutrons, Pp,n, are fixed by the conditions 

Z = J pp{v)dh , N = j PnWr . (18) 

Substituting the Hamiltonian density functional ^ in Eq. f|T7|) leads to the RTF-equations 
for protons and neutrons: 

+ ^7pP'° + + E*p^ = pp , (19) 
g^co' - gpP^' + E*p^ = p„ , (20) 



where Ep^^ = ^pj?^^ + rn*'^. For a spherical nucleus, the RTF equations fll9|20p together 
with the field equations fl8l)- f|TTl) completely determine the radial dependence of the proton 
and neutron densities and fields, i.e. pp^„(r), cr(r), u^{r), p^^{r) and ^°(r). The densities 

^ 1/3 

are obtained as the selfconsistent solution of Eqs. fl8l lTT|18ti20|) . with PFi = (3vr^pj) {i = 
p, n) according to the local density approximation. This method then yields a consistent 
description of the groundstate. The densities are realistic but somewhat too steep in the 
surface region. 

The solution of the RTF equations gives the nuclear density Pi{r) and the local Fermi- 
momentum PFii^) {i = P^n). The test particles for the numerical solution of the transport 
equation are distributed according to these functions. The propagation of the system is 
described by the test particle equations of motion, which directly follow from the transport 
equation ([2]) by putting its r.h.s. equal to zero and read (j = 1, ...,A ■ Af with M being the 
number of test particles per nucleon): 

= f (21) 

Pf = ^F'' + (22) 

At each time step of the simulation, the scalar density f|T2|) and currents f lT3|) - f|T5|) are 
calculated on a grid in coordinate space. Using these quantities, the equations of motion 
for the meson and Coulomb fields (I8i)-(|TT1) are solved numerically. Note, that the solution of 
Eq. dH]) requires some more iterations, since the scalar density ps in its r.h.s depends itself 
on the a meson field via the effective mass m*. 

The collision term in the r.h.s. of Eq.(|2]) is simulated by explicit two-body collisions 
between test particles using the geometrical collision criterium (c.f. 26|] for details). An 
important feature here is the numerical treatment of the Pauli blocking, which is discussed 
in detail below. 



B. Pauli blocking 

The frequency of two-body collisions in a Fermi gas depends on the occupancies of the 
scattering final states via the hole distribution functions fi (see Eq.(l6])). By the energy 
and momentum conservation, no collisions take place at zero temperature, when /j(r, p*) = 
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^(pi?- (r) — |p*|). In test-particle simulations, however, it is impossible to model the exact T = 
Fermi distribution. This causes some spurious two-body collisions even in the ground state 
nucleus. The magnitude of this spurious effect crucially depends on a numerical technique 
of the Pauli blocking calculation. In the standard GiBUU 29|, |38|, the occupation number 
/j(r,p*) is calculated by counting the number of test particles in the phase space volume 
element composed of small spherical volumes AVr with radius centered at r in coordinate 
space and AVp with radius Vp centered at p* in momentum space: 



where 



2AV..AVpM 4. /rf,,2exp|-^^ f24) 
(27r)3 (27ra2)3/2 j \ 2a^ ^ ^^^^ 





is a normalization factor. In Eq. fl23|) . the sum is taken over all test particles j of the type 
i = p,n whose momenta belong to the volume AVp. In the coordinate space, the test particles 
are represented by Gaussians of the width a cutted-off at the radial distance Vc- The default 
values of parameters are Vp = 80 MeV/c, = 1.86 fm, a = 1 fm, Tc = 2.2 fm. This set of 
parameters is a compromise between the quality of the Pauli blocking in the ground state 
and the smallness of statistical fluctuations in the case of simulations with A/" ~ 200 test 
particles per nucleon. Typically, this is good enough for modeling the heavy-ion collisions 
at the beam energies above ~ 100 MeV/nucleon. In the case of a small-amplitude dynamics 
near nuclear ground state, the accuracy provided by Eqs. fl23l) .( l24l) is not enough, when the 
default parameters are used. The main reason is the constant, i.e. momentum-independent 
radius Vp, which introduces a spurious temperature of the order of several MeV. To reduce 
this effect, we have introduced the coordinate- and momentum-dependent radius of the 
momentum space volume AVp as rp{r, |p*|) = max(20 MeV/c,pi7'.(r) — Ip"^]), which provides 
a sharper momentum dependence near Fermi momentum. The calculations of the present 
work are performed with Af = 1000 — 10000. This allows us to use the reduced parameters 
also in the coordinate space: = 0.9 — 1.86 fm, a = 0.5 fm, Vc = 1.1 fm. The calculations 
with the default and modified set of parameters for the Pauli blocking are compared in the 
following Sect. HVl 
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FIG. 1: (Color online) Time evolution of the binding energy per nucleon (panel on the top) and 
root mean square (rms) radius (panel on the bottom) for a ground state ^'''^Sn nucleus. Vlasov 
calculations using the (dashed) standard initialization and the (solid) improved initialization are 
shown. The filled circle in the top panel at t = fm/c gives the RTF- value of the binding energy. 



IV. RESULTS 



We study first the influence of the initialization method on the temporal evolution of 
nuclei in their ground states. The improved transport model is applied then to the dynamics 
of low energy nuclear excitations, which are simulated by initializing slightly expanded nuclei. 
We present results from pure Vlasov and full BUU calculations for different nuclei. If not 
indicated elsewhere, for the Vlasov calculations 10000 test particles per nucleon were used. 
The full BUU calculations were performed with 1000 test particles due to time limitations. 



A. Stability of the ground state 

Fig. [1] shows the time evolution of the binding energy per nucleon and the root mean 
square (rms) radius of a ground state ^"^Sn nucleus. The standard initialization method using 
the empirical Woods-Saxon density distribution produces the binding energy smaller by 0.3- 
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distance (fm) distance (fm) 

FIG. 2: (Color online) Density profiles of protons (panels on the top) and neutrons (panels on the 
bottom) for the same nucleus as in Fig. [TJ The thick curves are RTF calculations. The other 
curves show density distributions from the Vlasov dynamics at different times (as indicated) using 
the standard initialization (panels on the left) and the improved one (panels on the right). 



0.4 MeV/nucleon with respect to the RTF value of Eb/A ^ 8.1 MeV. This is expected, since 
the minimum of the total energy is not reached by the standard initialization. The binding 
energy varies with time due to numerical errors in the solution of the time evolution equations 
fl2T|) . f l22|) and field equations fl8|)-f lTT|) . For the standard initialization, the rms radius reveals 
quite strong fluctuations, comparable in the amplitude with the true GMR vibrations (see 
Fig. E]). These artificial temporal oscillations lead also to a significant particle loss with 
increasing time, if collisions are included (see below). Applying the improved initialization, 
in which the same Hamiltonian density functional is used for both the initialization of the 
nucleus and its temporal propagation, the situation becomes considerably better. At t = 
fm/c the value of the binding energy per nucleon agrees with the corresponding RTF value, 
and the rms radius stays almost constant in time. 

A more detailed picture of the Vlasov calculations with the standard and the improved 
initialization method is shown in Fig. |2] in terms of the proton and neutron density distri- 
butions. Using the standard initialization (figures on the left panel) the initial {t = fm/c) 
density profiles do not fit that one of RTF. This leads to significant density oscillations 
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FIG. 3: (Color online) The proton mean field potential + S (see Eqs. (OHI)) along the z-axis 
passing through the center of the ^"^"^Sn nucleus. The different curves show the Vlasov results at 
different times (as indicated) using the improved initialization (upper panel) and the standard one 
(lower panel). 



around the true ground state density profiles (RTF) with the result of spurious oscilla- 
tions in the rms-radius of the system (see again Fig. [1]). A consistent treatment between 
the ground state nucleus and its propagation leads to very good stable configurations, see 
graphs on the right panel of Fig. [2l 

We remind that in relativistic transport studies the central mean-field potential arises 
from the sum of the large negative Lorentz scalar and large positive Lorentz vector potentials. 
Thus, small spurious variations in density cause strong numerical fluctuations in the mean- 
field potential. This is demonstrated in Fig. [3l where the mean field potential is displayed 
as a function of the coordinate along the central z-axis. The Vlasov calculations with the 
standard initialization (panel on the bottom) show large fluctuations of the order of 10%, 
while these fluctuations almost vanish in the calculations using the improved initialization 
method. 

Including two-body collisions requires a careful implementation of the Pauli blocking to 
prevent the ground state to be destroyed. To give an impression of how well it is working 
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FIG. 4: (Color online) Momentum dependence of the proton occupation number /proton in the center 
of ^^C and ^'^'^Sn nuclei is shown in the upper and lower left panels, respectively. Radial dependence 
of /proton at the zero momentum for ^^C and ^'^''Sn is depicted in the upper and lower right panels, 
respectively. The results are presented for the standard (thin solid lines) and momentum-dependent 
(thick solid lines) radius rp. The fluctuations of /proton near Fermi momentum pp — 250 MeV/c 
are due to finite number of test particles per nucleon which was set to 10000 in this calculation. 
The nucleon density in units of po is shown additionally by dashed lines in the right panels. We 
see that at the half-central-density radius, the proton occupation number is only about 10% below 
unity. 



in our test particle calculations, Fig. S] shows the momentum (left) and radial (right) de- 
pendence of the proton occupation numbers, which are used in the evaluation of the Pauli 
blocking factors, for ^^C and ^"^^Sn nuclei. The calculation with the default parameters of a 
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FIG. 5: (Color online) The nuclear rms-radius Eq. (j26p (left panels) and the number of particles 
in the high-density region Eq. (j25p (right panels) as a function of time for the ground state ^^C 
and ^'^''Sn nuclei. The Vlasov calculations are represented by dashed lines. The full BUU results 
with the standard and modified Pauli blocking parameters are shown by thin and thick solid lines, 
respectively. 



Pauli blocking produces a rather diffuse momentum dependence, especially for the light ^^C 
nucleus. Using the momentum-dependent radius and the reduced width of a Gaussian, 
as explained in subsect. IIII Bt largely improves the momentum dependence of occupation 
numbers near Fermi momentum. The radial dependence of the occupation numbers also 
becomes closer to the step function, when the calculation is done with the modified Pauli 
blocking parameters. 

To demonstrate the effect of Pauli blocking parameters on the ground state evolution, 
in Fig. [5l we present the time dependence of the rms-radius of a nucleus (left) and of the 
number of particles in the high-density, p > pmin, space region (right) for the carbon and tin 
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nuclei. Explicitly, these quantities have been calculated as 



part 




(25) 



P>Pmin 



r-1 
part 




p(r) 



(26) 



P>Pmm 



with pmin = O.lpo- Vlasov calculations practically conserve the number of particles in the 
high density region and produce almost constant in time rms-radii. The spurious effect of 
two-body collisions in the ground state nuclei leads to the particle emission to vacuum which 
amounts after 200 fm/c to about 10% of the total mass number in the case of standard Pauli 
blocking parameters and ~ 5% in the case of the modified parameters. Correspondingly, the 
rms-radius gets reduced. This spurious reduction has to be excluded in the calculation of 
the rms-radius in full BUU simulations, as discussed later on. Below, we will always apply 
the modified Pauli blocking parameters as explained in subsect. IIIIBi 

B. Giant Monopole Resonance: GiBUU calculations 

The different methods of initialization of nuclear ground states influence the dynamical 
calculations of excited nuclei. We will model the low energy nuclear giant monopole collective 
excitations by initializing an expanded nucleus at t = fm/c. This is realized by re-scaling 
the coordinates of the test particles such that the corresponding excitation energy is close 
to the experimental values. The relation between the scaling parameter and the excitation 
energy is obtained by expanding the energy per nucleon E/A around saturation density or 
the ground state radius Rq: 



With = 258 MeV and using the experimental values ^] for the excitation energy 
= E — Eq one obtains scaling parameters AR/R = [R — Rq)/Rq in the range of 
~ 0.03 — 0.05 for A E (56,208). The system is then propagated either without (Vlasov 
mode) or with (BUU mode) collisions between the nucleons of the excited nucleus. 

Fig. shows the time evolution of a ^°°Sn nucleus using the different prescriptions of 
initialization. 




(27) 
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FIG. 6: (Color online) Time dependence of the rms-radius for the excited ^'^'^Sn nucleus. The 
meaning of the different curves is explained in the text. The calculations are done in the Vlasov 
mode. 



The standard method does not provide stable solutions and the rms-radius of the ^°°Sn- 
nucleus explodes for t > 200 fm/c due to strong particle loss. The situation is considerably 
improved in the new initialization. A clear oscillation signal can be seen with a particle emis- 
sion for very late times {t > 275 fm/c) which is, however, only moderate. Emitted particles 
increase the rms-radius of the total nuclear system which hinders the true oscillation signal. 
They have to be excluded, therefore, in the calculation of the rms-radius of an oscillating 
nucleus when extracting the excitation energy and width of the GMR. We thus consider 
only particles at densities higher than O.lpo according to Eqs. fl25|) . fl26|) . After this correction 
we obtain a clear signal for the resonance. The period of the oscillation characterizes the 
frequency and thus the excitation energy of the resonance, and an exponential damping (see 
below) — its finite lifetime. 

We have performed pure Vlasov mode calculations for different excited nuclei and ana- 
lyzed the results in terms of the time dependence of the rms-radius, as seen in Fig. [71 With 
increasing mass number the frequency of the oscillation decreases and thus also the exci- 
tation energy. An exponential damping is visible, even without the inclusion of coUisional 
effects. This effect has been interpreted as a wall friction 39, 4ol |. 



and it will be discussed 



later. An important feature in the study of nuclear collective excitations is the calculation of 
the lifetime or the width of the different multipole modes of the nuclear excitation. TDHF 
theory and Vlasov dynamics does not include any coUisional broadening effects. 
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FIG. 7: (Color online) Time dependence of the rms-radius for different nuclei, as indicated. (Solid) 
Vlasov calculations, (dashed) fit according to Eq. (p9|) . 



The nuclear collective dynamics of giant multipole vibrations in the ground state nuclei 
has not been so far investigated within the full BUU equation, mainly due to the reasons of 
a ground state instability. With the improved initialization method the Vlasov propagation 
of a ground state is almost perfect. The numerical procedure of the Pauli blocking in a full 
BUU ground state simulation is improved by applying the modifying method, but is still far 
from being exact. This situation leads to a spurious particle emission in ground state BUU 
calculations and thus to a spurious escape width Fe^, which impedes the determination of 
the total width. 

The spurious contribution to the total width has therefore to be excluded, as explained in 
Fig. [HI where the time evolution of bound particles for a nucleus in its ground state (dashed 
curve) and in the GMR mode (solid curve) is shown. Already in the ground state, particles 
leave the nucleus due to spurious two-body collisions resulting in a spurious escape width 
Fes. On the other hand, the BUU simulation for the excited nucleus contains also an escape 
width which is very similar to that of the ground state. Asumming Fes to be the same for 
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FIG. 8: (Color online) Time dependence of the number of particles in high-density region (top 
panel) and the r.m.s. radius (bottom panel) defined according to Eqs. (j25p . (|26p for a ^^Ni nucleus. 
Dashed lines: ground state BUU calculation, solid lines: BUU calculation for an excited ^^Ni 
nucleus in the GMR mode. 



both ground state and excited state we can make the following Ansatz: 
Ngs(t) = A^o exp(— Test) for the ground state 

Ne^^{t) = Noexp{-T^st)F{t) for the excited GMR state , (28) 

where Nq is the number of particles at t = fm/c and F{t) is a fit function for the oscillation 
signal, which contains the stochastic collisional width: 

F{t) = a + 13 cos(wt + 6) exp(-7t) . (29) 

The spurious width can be excluded by taking the ratio of Ngs and N^xc, or their difference in 
the case of very small value of Fes- The latter method is applied here in extracting the GMR 
width, since Fg^ ~ 0.044, 0.045 MeV for a ground state and an excited nucleus, respectively. 
Thus, the physical escape width (relative to that in the ground state) is almost negligible, 
and we obtain the total width in full BUU calculations by a fit according Eq. (!29l) to the 
so-called corrected r.m.s. radius defined as 

< rLr it) =< rL >^/^ it)- < r% >''^ (t)+ < rj, >^/^ (0) , (30) 
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FIG. 9: (Color online) The corrected r.m.s. radius of Eq. ()30p plotted as a function of time for 
different nuclei — solid lines. Calculations are done in the full BUU mode. The fit according to 
Eq. ([29|) is shown by the dashed lines. 



where the depletion due to spurious particle loss is subtracted. However, the damping rate 7 
of the corrected r.m.s. radius oscillations is still somewhat influenced by an imperfect Pauli 
blocking procedure, which makes some uncertainty in our calculations of the GMR width. 
We also have checked the consistency of the extracted results on E* and F obtained with 
the subtraction method by performing additional Fourier analyses of < >^/^ (t) before 
and after the correction. Lorentzian fits to the Fourier spectra lead to values for the width, 
which are the same with those extracted with the subtraction method. Also the excitation 
energy is not essentially affected. However, the excitation energy in full BUU calculations 
is slightly above the corresponding values in Vlasov mode. This effect is due to the spurious 
particle emission in full BUU calculations, which is furthermore related to the collective 
response of a smaller system. 

The results for the corrected r.m.s. radius produced by full BUU calculations for different 
nuclei are shown in Fig. |9l Again, a clear oscillation signal is visible, however, due to the 
inclusion of collisions, a significant damping of the oscillations of the rms-radii appears, as 
compared to the pure Vlasov calculations of Fig. [71 We note, that in the case of a pure 



20 



Vlasov mode calculation, the difference between corrected and not corrected r.m.s. radii is 
negligibly small. 

We have performed a fit to the results displayed in Figs. [7] and M according to Eq. ( 12^ 
to obtain the centroid energies and the total width of the breathing mode. The results are 
summarized in Fig. [TO] in terms of the excitation energy E* = u and the damping width 
F = 27 of the GMR, as a function of the mass number. The factor of two in the last 
formula is motivated by our intention to report the Full Width at Half Maximum (FWHM) 
of the GMR strength 17|. The gray band in the transport calculations for the width is 
an estimation related to the numerical uncertainty of the Pauli blocking factors. We first 
discuss the results obtained for E* and then those for the width. 



C. Giant Monopole Resonance: discussion 

The theoretical Vlasov-calculations for the excitation energy scale with A~^^^ and over- 
estimate moderately the experimental data. The mass dependence is also consistent with 
TDHF calculations and with microscopic RPA studies Similar results for the excitation 
energy are obtained using full BUU calculations. It is a well known question, why relativistic 
structure calculations can explain the excitation energies of GMR with a higher value of the 
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42I) tha 



incompressibility modulus K^q = 250 — 270 MeV (c.f. refs. 
MeV deduced from nonrelativistic approaches It has been shown in Ref. 



han ^ 220 
4l[ |. that this is 



at least partly related to a stiffer density behaviour of the symmetry energy in RMF models 
with resp ect to the Skyrme-type effective interactions. On the other hand, the calculations of 
Ref. [lOj have revealed the influence of differences in the surface compressibility in different 
RMF models on the GMR frequency. There is also another possible reason for differences 
between RMF and nonrelativistic approaches, which we address below. 

Our calculations have been performed without nuclear Lorentz force, i.e., without taking 
into account the space-like components of the vector field. 

We expect that a nonrelativistic Vlasov calculation, employing the same energy-density 
functional suitably parameterized, e.g. in a Skyrme form, will produce very similar results. 
To study the influence of the nuclear Lorentz force on the excitation energy (and the width) 
of the GMR, a Vlasov calculation explicitly taking into account the space-like components 
of the vector field has been performed for the lead-208 nucleus. The results are shown 
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TABLE I: Excitation energy E* for a ■^''^Pb nucleus in different models: (Vlasov wLF) Vlasov 
calculation including the Lorentz force, (Vlasov) Vlasov calculation without the Lorentz force, 
(RPA) RPA calculations by Lalazissis et al. |32]. All the results are given in the case of NL3* 
model (i^oo = 258 MeV). 



model 


Vlasov 


Vlasov wLF 


RPA 


exp. 


E*" (MeV) 


15.2 


14.2 


13.9 


13.7± 0.5 



in table U and compared with other theoretical models of nuclear structure. First of all, 
the effect of the Lorentz force is negligible for the GMR width (not shown in the table). 
The excitation energy is moderately affected. In particular, a decrease of E* by ~ 6.5% 
towards the experimental data and the relativistic structure calculations is observed, when 
the Lorentz force is included in the Vlasov calculation. We note that the same value for the 
compression modulus has been used in both Vlasov calculations. This result indicates the 
importance of the genuine relativistic effects when extracting the incompressibility modulus 
from GMR studies. 

Another feature of interest in the transport results using both modes, full BUU and 
Vlasov, is the moderate decrease of the monopole frequencies with increasing neutron excess, 
as one can see in the results for E* of Fig. (upper panel) for the Sn-isotopes. Such a 
trend is also supported experimentally |43]. It is well known that the monopole frequencies 
are affected by the slope of the symmetry energy around saturation, as in detailed discussed 
in Ref. It would be a challenge to extend this transport study to particularly isospin 

asymmetric systems, such as the Sn-isotopes, where systematic experimental studies exist 



431], and to investigate more exotic collective modes in neutron-rich systems. 

To understand the GiBUU results on the mass dependence of the GMR centroid energy 
on a qualitative level, we have performed some simple estimations. As it is well known from 
empirical GMR systematics pj, the centroid energy of the GMR follows the A'^/^ law: 

E* = 7]A-^/^ . (31) 

This behaviour can be understood as a consequence of a sound-like excitation in a finite 
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system, i.e. E* = Vsk, where Vg is a sound velocity and k = n/R is the eigenvalue of 
the lowest compressional mode determined from the disappearance of the pressure at the 
free surface {45], R = 1.2A^^^ is the nuclear radius. The hydrodynamical model 45] would 
give Vg = (Kcc/9mY^'^ ~ 0.17, which is the first sound velocity. Here, we have used the 

value of the incompressibility modulus K^q = 258 MeV provided by the NL3* model. The 

1/2 

dependence rj oc K/ , where Ka is, however, the incompressibility modulus of the finite 
nucleus, is also provided by the scaling model of the GMR . This model is usually applied 
in extraction of from experimental data on GMR energy (c.f. H) by using the relation 



E* = a/ Ka/tti < >, where < >~ 3R^/5 is the r.m.s. radius of the nucleus. 

On the other hand, according to the Fermi liquid theory [46|] the low-temperature col- 
lective excitations in the infinite system are of the zero sound type. It is well known, that 
the propagating zero sound type solutions of the dispersion relation for coUisionless Fermi 
liquid at zero temperature exist only for the repulsive particle- hole interactions (c.f. 
which is not the case for the NL3* interaction used in the present work (see below). Intro- 
duction of finite temperature, generally, restore the collectivity for attractive particle-hole 
interactions, although it was proved only for momentum- independent interactions [47]. It 
is quite difficult, therefore, to actually identify the GMR vibration as the zero sound mode. 
Assuming, nevertheless, the nuclear matter zero sound nature of the GMR vibration, we 
can estimate the sound velocity as ~ = pp/m*^ ^ 0.42, where pp = 257 MeV/c is the 
Fermi momentum at the nuclear matter saturation density and = \/ m*"^ + = 0.65m 
is the Landau effective mass in the case of the NL3* model. 

As a result, one gets the values 90, 111 and 223 for the coefficient r] in Eq. (|3Ti) for the 
hydrodynamical, scaling and zero sound pictures of the GMR, respectively. In the case of the 
scaling model, we assumed that Ka = K^o, which is a quite rough assumption. Generally, 
one has Ka < K^ mostly due to the surface contribution 1^. It turns out (see the upper 
panel of Fig. [TOl) . that Eq. ( !3TI) with the "hydrodynamical" value of = 90 well fits the 
Vlasov results for large masses numbers A > 100, although the reason for this is not fully 
clear for us. 

For light nuclei, as one can see from Fig. [TOl the transport calculations overestimate the 
experimental data on E*. This is the region where the RTF method becomes unreliable 
because the surface properties are not well described. However, the experimental determi- 
nation of the GMR parameters in light nuclei is rather uncertain due to strong fragmentation 
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FIG. 10: (Color online) Excitation energy (upper panel) and width (lower panel) of the GMR as 
a function of the mass number. Vlasov and full GiBUU results (as indicated) are compared with 
experimental data (open diamonds). The solid line in the upper panel shows the fit by Eq. (j3ip 
with 7] = 90. In the lower panel, the dashed line depicts the fit of the Vlasov calculation using the 
wall formula ()33p with = 20, while the solid line shows the full width F = —2a;/ (see Eq. ()35p ) 
taking into account both the wall and collisional contributions. Experimental data are from Ref. 
48|. 



of the 0^ strength 

The situation for the width is more involved. Pure Vlasov calculations predict a very 
small value for the GMR width and do not fit the data, as expected. The inclusion of 
collisions in the full BUU calculations improves the comparison between theory and experi- 
ment considerably. The underprediction of the theoretical calculations to the data becomes 
smaller, but an exact agreement is not achieved. For a deeper interpretation of the GiBUU 
results on the mass dependence of the GMR width several additional analytical calculations 
were performed. 

The width of a collective vibration is related to the dissipation processes. Within a pure 
mean field Vlasov calculation, the only damping mechanism is the one-body dissipation 
governed by a coupling of the single-particle and collective motions. Specifically for finite 
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systems, damping arises due to collisions of particles with a moving wall. This one-body 
dissipation leads to the following "wall" formula for the collective energy dissipation rate 



39| 



E = mpv j) n^da , (32) 

where p = pp + pn is the nucleon density, n is the normal component of the wall velocity, and 
the integral is taken over the surface of a vessel, v is the average speed of particles in the 
vessel. For cold nuclear matter one has v = ^vp. The formula (15^ is derived under strong 
simplifying assumptions of the gas at rest inside the vessel and of a sharp potential wall 
driven through the vessel. These conditions are usually assumed to be valid for small surface 



vibrations of the incompressible nuclear droplet j39|, possibly with some modifications [40 1 
(see also [l7| and refs. therein). Application of the wall formula to the compressional modes 
is more questionable. Nevertheless by performing a rather simple calculation for the GMR 



mode within the liquid-drop mode 
the one-body relaxation time (c.f. 



with a free surface we arrived to the usual formula for 

m) 

Twall = , (33) 

V 







with ^ = 0.5 which means an extremely strong dissipation. E.g., for ^°*Pb we obtain the 
wall dissipation contribution to the GMR width r^aii = 2/r„aii — 17 MeV. This result has 
to be considered as an extremely rough approximation. The self-consistency corrections, i.e. 



the collective motion of nucleons near the surface region |40| will modify Eq. fj32|) and r„aii- 
Other effects like the surface diffuseness and curvature are not taken into account by the 
wall formula at all. Having these reservations in mind we will treat ^ in Eq. f l33|) as a free 
parameter and determine it from comparison with the results of calculations in a Vlasov 
mode. This produces the value ^ ~ 20 as demonstrated in the lower panel of Fig. [THl 

We will now discuss the two-body dissipation. It follows from the Uehling-Uhlenbeck 
collision integral that the two-body collisions take place only in the case of local deviations of 
the Fermi surface from a spherical shape or/and in the case of a finite temperature. Moreover, 
the relaxation rate of a nonequilibrated Fermi gas toward thermal equilibrium depends 
mainly on the total excitation energy and not on the concrete shape of a Fermi surface 
deformation j49|. In particular, in the linear approximation with respect to the deviation of 
a distribution function from the local equilibrium the relaxation rate is proportional to 
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50| . This means, that coUisional damping of a small Fermi surface distortions practically 
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vanishes at T = 0. 

Thus colhsional damping of a vibrational motion excited in the ground state nuclear 
system is practically absent at the beginning of time evolution and gradually switches on as 
some part of the collective vibrational energy is transferred to the heat. In Refs. [l^ ll^ . 
an "apparent temperature" has been introduced in the calculation of the collisional widths 
of the giant quadrupole and giant dipole vibrations built on the ground state nuclei. The 
"apparent temperature" has been extracted in Is], Q by subtracting the collective energy 
-EcoU of a vibrational mode from the total excitation energy E*: T = a/ {E* — Eco\\)/a, where 
a is a level density parameter. This "temperature" has the real physical meaning only when 
the system reached a complete thermal equilibrium. Since the Pauli blocking factors in the 
Uehling-Uhlenbeck collision integral depend crucially on a temperature, calculations taking 
into account the time- dependent "apparent temperature" strongly increase the spreading 
width of a collective mode jisl . \\\ to a good agreement with experiment. In our estimates 
of the two-body dissipation we will use the upper limit for the "temperature" by putting 
Ecow = 0. For the level density parameter we will use the Fermi gas expression a = tc'^A/AEf. 

The GMR mode can be considered as a sound-like excitation inducing the Fermi surface 



2 



distortions of all multipolarities I > 2 [17[. At T ^ Ep the relaxation time of a small 

3) ri = K,/r 



amplitude Fermi surface distortion with multipolarity / > 2 is (c.f. Q, 
where ki depends on the NN cross sections. The parameters ni have been computed in 19l.l51| 
both for isovector and isoscalar vibrations for various choices of the NN cross sections. In the 
case of realistic energy and angular dependent vacuum NN cross sections this resulted in the 
following values of the isoscalar parameters: ki = 868, 881 and 401 MeV^fm/c for / = 2, 3 
and -|-oo, respectively. The case / = -|-oo corresponds to the relaxation of a single particle- 
hole configuration. For simplicity, we will set the same collisional relaxation time Tcou for all 
multipolarities / > 2: Tcqu = k/T^, where k = 3/(^2"^ + + k^) = 628 MeV^fm/c. In the 
limit of large relaxation times the interrelationship between various sources of dissipation 
can be ignored and one can treat also the wall dissipation as an additional source in the 
collision term by defining the total relaxation time as 



T 



- ' - r-oi + r^aii • (34) 



Here we neglected the particle emission from an excited nucleus. 

Applying the formalism of the linearized Landau- Vlasov equation in the relaxation time 
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approximation 17|, |25| leads to the following approximate expression for the imaginary part 



uji of the giant multipole resonance frequency u = ujr + iuj : 

uji ~ -qujR — — , (35) 

where the factor q is related with the Landau parameter Fq a.s q = 2/5(1 + Fq). Eq. fl35]) 
is, in fact, a suitable interpolation between the two well known limits |46|] of rare collisions 
Urt ^ 1 (zero sound) with uj ~ — 1/t and of frequent collisions Urt ^ 1 (first sound) with 
ujj ~ -qujRT. 

The Landau parameter Fq can be expressed via the nuclear matter incompressibility Koo 
and Fermi energy Ep = p|./2m2 as K^o = QEf{1 + Fq). The Landau effective mass is 
connected to the Landau parameter Fi as = m (l + Using the NL3* parameter 

set of the RMF model we obtain Fq = -0.20 and Fi = -1.04. 

We have applied Eq. ( I35l) with ur = E* to compute the GMR width F = —2ujj. The 
result is shown by the solid line in the lower panel of Fig. [TOl We observe that Eq. ( 135|1 
gives a reasonable estimate of the magnitude of the collisional broadening for heavy nuclei. 
However the GMR width computed by using Eq. (!35|) is too large for medium and light 
nuclei. This is expected, given the fact that we have used the upper limit of the "apparent 
temperature" , which becomes unphysically high for small mass numbers. 

We have to point out that the zero sound damping conditions are valid for medium- 
to-heavy mass region A > 50, where we have urt > 6. This creates a puzzle, since, 
as we have seen above, the nuclear matter zero sound model would strongly overestimate 
the experimental GMR centroid energies. The answer could be, that the finite size effects 
essentially modify the zero sound mode in a real nucleus. 

The GMR damping mechanism has been a long-winded problem in quantal structure 
calculations in the spirit of the RPA [52!] and GCM calculations 10|. It is not the scope 
of this work to list all the various structure calculations and discuss their details, however, 
they can serve for a qualitative comparison with our calculations. For more details we refer 
to review articles 53|, |5J]. Fig. [11] displays again the GMR width in the Vlasov approach 



and the full BUU transport model and shows a comparison with the two sets of calculations 

n 

from [55!]. 

The (Q)RPA calculations, except for the contribution from pairing effects in open-shell 
nuclei, can be considered as the quantum analogue of our semiclassical Vlasov-mode cal- 
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FIG. 11: (Color online) Width of the GMR as a function of the mass number. Vlasov and full 
GiBUU results (filled circles and gray band, respectively) are compared with experimental data 
(open diamonds) taken from [48]. The solid-circle and dashed-square curves shows the (Q)RPA 



and QTBA result, respectively, taken from Ref. 



55|. 



culations. Collective response is a superposition of Ip-lh excitations in the both types of 
calculations. The difference between (Q)RPA and Vlasov results for the GMR width stems 
from missing the quantum fragmentation (Landau damping) width contribution in our cal- 
culations^. Including a quasiparticle-phonon coupling in the quasiparticle time blocking 



56|, 



approximation (QTBA) |55| can be regarded as a coupling to the 2p-2h configurations 
which leads to the strong increase of the total GMR width. We observe a similar effect 
in our full BUU calculations, since nucleon-nucleon collisions in the Uehling-Uhlenbeck col- 
lision term also generate 2p-2h excited states. Moreover, the difference between QTBA 
and (Q)RPA results is quite close to the difference between full BUU and Vlasov results. 
This again indicates the importance of the correct description of the fragmentation width 
contribution. 

In the spirit of refs. jisl. IigIi. where the authors argue, that the Markovian approximation 
(i.e. a standard Uehling-Uhlenbeck collision term) is not able to produce any broadening of 
the giant multipole vibrations at zero temperature due to severe restrictions of the available 
phase space for two-body collisions, our result on the width enhancement by two-body 



^ In the analytical models based on BUU equation [l3, [l3 the quantum fragmentation width contribution 
has been included as an additional source term in kinetic equation. However, this is a pure phenomeno- 
logical way to describe the damping width of giant multipole resonances in the ground state nuclei. 
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collisions looks quite surprising. In [16|], the authors explain about 25-30% of the observed 
GMR width by taking into account the memory effects in collision term (non-Markovian 
approach). However, the analytical models ISMlTl. Il9| are based on the linearized kinetic 



equation violating the energy conservation. In particular, the temperature increase due to 
the damping of collective motion is completely neglected in these approaches. We stress, 
that the solution of the full non-linear BUU equation in finite system leads to much stronger 
collisional broadening, than expected from the linearized models. 



V. CONCLUSIONS AND OUTLOOK 



The aim of the present work has been to study theoretically the excitation energy and 
the width of the giant monopole resonance in the framework of a semiclassical transport 
model. For this purpose, a description of nuclear ground states and their phase-space evo- 
lution within a unified semiclassical framework is indispensable. Nuclear ground states are 
described by the RTF method using a relativistic Hamiltonian density functional which, in 
particular, contains the space-derivatives of the mean meson fields. The proton and neutron 
densities from the RTF calculations serve to generate initial test particle configurations of 
different nuclei. The time evolution of the nuclear system is calculated by solving the kinetic 
equations which employ the mean-field potentials consistent with those used in RTF. The 
improved transport model gives a perfect stability of different ground state nuclei over long 
time scales in pure Vlasov dynamics. 

The improved initialization method was found to be important to generate a clear signal 
of the breathing mode. Except for light nuclei, pure Vlasov calculations predict a mass 
dependence of the excitation energy which is consistent with available experimental data and 
with a simple liquid-drop model. The situation is, however, more complex concerning the 
lifetime of the breathing mode. Vlasov simulations strongly underpredict the experimental 
data on the GMR spreading width. The GMR width calculated by using the Vlasov equation 
behaves as cx A^^^^, which is consistent with the (modified) wall formula. 

In order to better understand the GMR damping mechanism, the full transport calcula- 
tions of the GMR mode have been performed for the first time. The inclusion of two-body 
collisions strongly enhances the total GMR width to a better agreement with experimental 
data. The strong damping of the GMR including the collision term can be understood in 
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terms of analytical models for one-body and two-body dissipation taking into account the 
temperature increase due to the dissipation of the GMR motion. 

The Pauli blocking strongly influences the dynamics in full BUU calculations. Thus it 
has been treated as precise as possible in the present work. The spurious particle emission 
due to incomplete (numerical) Pauli blocking destroys the stability of nuclei on the long 
time scales of the order of ten periods of GMR oscillations and has been subtracted. This, 
however, produce a systematic error of about 30% in our results on the GMR damping 
width. 

Overall, the full BUU calculation underestimates the total GMR width by about 30—50%. 
This might be related to the missed Landau damping contribution in our scmiclassical 
approach, as the comparison between Vlasov and RPA calculations may indicate. Moreover, 
the memory effects in a collision integral neglected in our calculations also increase the widths 
of giant multipole resonances. 

In conclusion, the breathing mode energies and widths in medium-to-heavy nuclei are 
reasonably well described within the GiBUU approach, when the same Hamiltonian en- 
ergy functional is consistently used in the initialization procedure of nuclear ground states 
and their phase-space evolution. Thus, an extension of the present work to investigate 
other modes of collective excitation, such as isovector dipole and isoscalar quadrupole res- 
onances, seems possible and would be a helpful tool to understand better the dynamics in 
low energy reaction physics. Future transport applications to fusion/deep-inelastic collisions 
for isospin asymmetric systems to investigate exotic collective modes in neutron-rich finite 
systems, such as the pigmy dipole resonance, seem possible. Furthermore, the improved 
transport model may be a better tool in theoretically describing hadron-induced reactions, 
in-particular the proton-induced those leading to nuclear fragmentation, and very peripheral 
heavy-ion collisions, in which the stability of nuclear ground states is essential. 
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